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Abstract. The main implications of noncommutativity over astrophysical objects are 
examined. Noncommutativity is introduced through a deformed dispersion relation E 2 — 
p 2 c 2 (1+A-E) 2 +m 2 c 4 and the relevant thermodynamical quantities are calculated using the grand 
canonical ensemble formalism. These results are applied to simple physical models describing 
main-sequence stars, white-dwarfs and neutron stars. The stability of main-sequence stars and 
white dwarfs is discussed. 



1. Introduction 

Noncommutativity is believed to be an important feature of space-time at quantum gravity 
scales [Tj[2]. Interestingly, it also appears in the context of string theory [31 H]. 

Noncommutative quantum field theory (NCQFT) models can be implemented by substituting 
the normal product between fields by a noncommutative (NC) one, the so-called Moyal 
product [5j [6]. This procedure introduces a minimum length scale, but does not solve, as 
one could expect, the renormalization problem of the quantum field theory. One also encounters 
additional difficulties on issues such as causality and unitarity. Related questions associated 
with noncommutativity involve the violation of translational invariance [7J, scalar fields and 
their stability in curved spaces [8l [9J [10] . NCQFT can also be implemented by generalizing the 
Heisenberg-Weyl algebra of noncommutative quantum mechanics [HI [121 E21 EH ESI El EZ1 EE] 
to field algebra [TTj. This approach can be used, for instance, to investigate the inflationary 
period [19]. Another way to introduce noncommutativity in field theory is through a deformed 
dispersion relation in a commutative space-time rather than considering a noncommutative 
space-time over which the fields are defined [20]. Since there is no consistent theory of 
noncommutative gravity |21|, [22] this approach is justified. This is inspired by quantum groups 
methods [23J, by a possible breaking of Lorentz symmetry at high-energies |24[ I25 [ [26] and the 
modifications of special relativity that take into account the possible existence of a minimum 
length [271 [28]. 

One mentions here some possible ways on how noncommutativity can lead to a deformed 
dispersion relation. The first approach is related to quantum groups methods [23]. Assuming 
the following commutation relations 



[x\t]=i\x i [x\x j ] = 0, (1) 

where i,j = 1,2,3 and A is constant, it can be shown that the associated relationship between 
dynamical quantities is given by [23] 

A- V S + e~ XE - 2) - p 2 c 2 e- XE = m 2 c\ (2) 

which is a deformed dispersion relation. 

A relationship between noncommutativity and a deformed dispersion relation also arises in 
the context of special relativity proposal with an invariant length |27j . In this instance, one 
modifies the boost generator and from the invariant quantity 

where Lp is the invariant length, a deformed dispersion relation follows. The thermodynamics 
of an ideal fluid following the deformed dispersion relation Eq. ([3|) was investigated in Ref. |29j . 

From the previous discussion, it is natural to expect that NC arises at high energies, and 
hence in the early universe [20\ 130] . However, in this contribution one discusses the results of 
Ref. [31] where NC low-energy effects on astrophysical objects are investigated. Also here the 
noncommutativity is introduced via a deformed dispersion relation. Using the grand canonical 
ensemble, one finds the leading noncommutative correction to particle number, energy densities 
and pressure. This method is applied to radiation, non-relativistic ideal gas and degenerate 
fermion gas and then used to describe the main features of main-sequence stars, white dwarfs 
and neutron stars. 

2. Deformed Dispersion Relation 

One can define a simplified version of Eq. ([5]) as in Ref. [2D] 

E 2 = p 2 c 2 (l + \E) 2 +m 2 c\ (4) 

Indeed, performing A — > 2A in Eq. ([5]), Eq. ^ matches Eq. up to the first order in A |2U| . 
Solving for E one gets 

Xp 2 c 2 + a/p 2 c 2 + m 2 c 4 (l — A 2 p 2 c 2 ) 
E = 1 - AVc 2 ' (5) 

The general behaviour of Eq. ([5]) is depicted in Fig. [TJ One considers the particle 
branch, where all energies can be attained although there is a maximum momentum given 
by Pmax = 1/ Ac. Thus the parameter A is associated with the maximum momentum |3U] . Since 
one aims to investigate energy scales encountered in astrophysical objects and A is presumably 
related to the inverse of the quantum gravity energy scale, it suffices to keep only the first 
correction in A. Therefore, at first order in A, Eq. ([5]) reads 



E = \p 2 c 2 + yp 2 c 2 + m 2 c 4 . (6) 
Notice that the usual relativistic dispersion relation is recovered for A — > 0. 
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Figure 1. Deformed dispersion relation for the particle branch of Eq. ([5]). 



3. Deformed Statistical Mechanics 

The foundations of statistical mechanics are not modified for more general dispersion relations 
|20t 1321 [33] . One considers the grand canonical ensemble, since it is the suitable formalism to 
describe radiation, ideal gas and degenerate fermion gas. Consider a system with ./V particles. 
Each state is labeled by j (j = 1, 2, • • •) and corresponds to rij particles with energy Ej. One 
defines the fugacity as e^, where /j, is the chemical potential, (3 = (fcgT) -1 , ks is the Boltzmann 
constant and T the absolute temperature. The grand canonical partition function is defined as 

2=eii [*r pE >] nj = e n [^ Ej) ] nj ■ ( 7 ) 

n 3 3 n J 3 

Through the grand canonical potential, one obtains the thermodynamics 



$ = -PV = --\nZ = — l - M 1 + aze-^) (8) 
P a P j 

where P is the pressure, V the volume, a = 1 for fermions and a = — 1 for bosons. Using Eq. 
Q, the pressure is obtained in the large-volume limit X^_e — ^ I (2nK$ > 

where E = E(x,p) and 7 is the multiplicity of states due to spin. At this point one has to 
introduce an specific relationship for E = E(x,p) and one chooses Eq. ([5]) for this function. Since 
the deformed dispersion relation Eq. ([5]) depends only on the absolute value of the momentum, 
one can compute the configuration variable integral in Eq. ([9|), which yields the volume V. The 
pressure then reads 



2vr 2 ^ 3 Jo \zr 1 eP E + a) \3 dp J 

Notice that the upper limit of this integral is p max = (Ac) -1 . The average number of particles 
and the average energy are given respectively by 

^ = -(tt) =E»( js 0. <JS>=5>n(JE0, (11) 

where n(E) = (^z~ l e^ E + is the occupation number of particles with energy E. Considering 
the large-volume limit and integrating, one finds the particle number and energy densities 
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d P -I^B i - > ( 12 ) 



V 2vr 2 ?i 3 7o I z- x eP E + 
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U= V = ^ = ^t dPp2 (z-Jz + a)- (13) 
3.1. Deformed Radiation 

A gas of photons is used to describe radiation. For photons, the deformed dispersion relation 
Eq. (gj) reads 

E = pc(l + XE). (14) 

To obtain the energy density one substitutes Eq. (TT4"|) into Eq. (|13p with a = — 1, fi = and 
7 = 2, that is: 

1 [°° E 3 I 
U ~^¥^Jo dE e^-l(l + XEr (15) 

This integral cannot be solved analytically. The approximation XksT C 1 is introduced. 
This limit is justified for astrophysical objects since the highest temperatures attained in stars 
are around (10 11 — 10 12 ) K [34J, yielding the condition A <C 10 GeV -1 which is satisfied if 
A ~ Eqg w here Eqq is the quantum gravity energy scale. The limit, Xk B T 3> 1, presumably 
relevant at the early universe, was investigated in Ref. |20j . 

Defining a variable y = XE. Eq. (|15p can be written as 

1 f°° y 3 1 

dy — — — TT— -u- (16) 



vr 2 ft 3 c 3 A 4 Jo y e T^r _ 1 (l + y) 

As XksT <C 1, it suffices to expand y 3 /(l + y) 4 in Taylor series around zero, since away from 
the origin other values will be exponentially suppressed by (e Afc s T — One getdU 



1 
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where a 



7r 2 ft 3 c 3 A 4 

7T 2 fc 4 
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96C(5) 
ir 2 h 3 c 3 



Xk%T 5 



(17) 



60h s c 2 



is the Stefan-Boltzmann constant. Notice that the Stefan-Boltzmann law 



is recovered for A — > 0. This expression can be rewritten as u = 4o" e g(T; A)T 4 /c, where the 
effective Stefan-Boltzmann "constant" is given by: 



i I 1440C(5) 
o" eff (T; A) = a I 1 ^ Xk B T 



7T> 



(18) 



So the leading noncommutative correction reduces the energy density. Using Eq. (TIT)]) , the 
pressure is 



1 



3vr 2 n 3 c 3 Jo 



oo £3 
dE 



1 



- 1 (1 + AE 1 ) 5 



(19) 



This integral can be solved in the same way as the energy density integral. One obtains 

1 Using the formula J °° ' eM , J 1 * = -4 7 r(i/)C(i / ) [Re/i > 0, Re;/ > 0], where T(u) and are gamma and zeta 
functions, respectively [35] . 
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•J JO e \k B T _ l JO e Xk B T "JC 7T z n C J 



4<r 4 24C(5) 



7T 2 n 3 C 3 A 4 

To find an Equation of State (EoS), one must divide Eq. (|20|) by Eq. (|17|) 



+ (Xk B T)^p-, (21) 



£ = I + (Ater) ^) 
u 3 7r 

Using Eq. (fT7|) . one writes T as a power series in A. Substituting this result into Eq. ([21 
one finally gets the leading noncommutative correction to the EoS 



P = ^ 
3 



1 360C(5) fl5\ 1/4 



Notice that the usual relationship u = 3P is recovered for A — > 0. 



(22) 



3.2. Non-relativistic Ideal Gas 

Consider Eq. ([6]), hence for a non-relativistic ideal gas p -C mc, it reads 

E = Xp 2 c 2 + mc 2 + ^ + o(^-)\ (23) 
2m \mc J 

The condition to be satisfied in order to ensure the relevance of the noncommutative correction 
is A ^ (2mc 2 ) -1 . However, for matter found in main-sequence stars, hydrogen for simplicity, 
one finds A < 1 GeV -1 , which certainly is a too stringent bound. Hence, one does not consider 
any noncommutative correction to the non-relativistic ideal gas. The pressure is the well know 
expression [32] P = NksT/V . After introducing the mean molecular weight /j>n = p/(nm^), 
where n = N/V, p is the mass density and m^ is the nucleon mass^l, the pressure reads 

P = — - — k B T. (24) 

3.3. Deformed Degenerate Fermion Gas 

A fermion gas becomes degenerate at low temperatures (T — > 0). In this case, the fermionic 
character gives origin to a very simple expression for the occupation number. In terms of the 
momentum, the occupation number is n(p) = H(pp — p), where pf is the Fermi momentum, 
above which there are no occupied states and H{x) is the Heaviside step function. For stars, 
the momentum varies between MeV/c — GeV/c, so pf <C (Ac) -1 and one uses Eq. (0) as the 
deformed dispersion relation. 

For spin one-half particles (7 = 2), the particle number density Eq. (|12|) is given by 

N 1 fT-c 2 1 [pf 2 p 3 F (mcfx 3 

n= v = ^Jo dpn{p)p= ^Jo dpp= ^K 3= ^T' (25) 

where x = pF/mc. Using Eqs. (0 and (fTHj) . one obtains the energy density 

„2 




p 2 c 2 + m 2 c 4 H — —^7 / dp p 4 

7T 2 ?T Jo 

M + (Amc 2 )^) . (26) 



2 m N = lamu = 931, 494 MeV/c 2 [36] 



where f(x) = x{2x 2 — 3)vl + x 2 + 3sinh 1 x [37]. This shows that the first noncommutative 
correction is given by Xmc 2 . The pressure is obtained using Eqs. ([6]) and (|10|) : 



c 2 /"^ , p 4 2\c 2 [Pf , 4 (mc 2 ) 4 //(a;) 2 ^ 5 \ , „x 

4. Application to Astrophysical Objects 

^.i. Main-sequence Stars: The Sun 

The simplest model for a main-sequence star consists in assuming that the star is composed 
by a non-relativistic ideal gas and radiation that are maintained in equilibrium by gravity |37j . 
The total pressure P comprises the sum of the radiation pressure P rac i (Eq. ([20)) ) and the gas 
pressure -Pgas (Eq. ((21]) )■ One denotes the ratio between Pgas and P by Ps- One then has 
P = Pgas /As = f ra( j/(l — Ps)- Substituting Eqs. ([20]) and (l23|) into this relationship one finds 

" -*bT=^^-T±-X J:!^ Wl \T\ (28) 



Hn™nPs 3c(l-/3 5 ) vr 2 ?i 3 c 3 (l-/3s)~ 

The aim is to obtain P = P(p), so one must solve this equation for T, obtaining an expression 
T = T(p) up to the first order in A. Inserting this result into Pgas = /3P and using Eq. ([21 
one finds that P = ifi/9 4 / 3 + AK 2 p 5 / 3 where 



1/3 



(l-fe) 1 / 3 / 3cfc 4 , 
(/3s^v) 4 / 3 \*rm$J 

(1-&?) 2 / 3 / 360C(5) \ 
((3sPn) 5/3 V vr 4 J 



2.67 x 10 



10 



2/3 



4am 



5/2 

AT ' 



(1-/35) 1/3 
(/%^) 4 / 3 



4.52 x 10" 6 



J m 



kg 



4/3 



(1 ~ fe) 2/3 
(/%Wv) 5 / 3 



J 2 m 2 



5/3 



kg 



(29) 



(30) 



One assumes the Eddington's Standard Model of Stars [33]: (3s is constant inside the star, 
and pj\[ is constant since the chemical composition is not supposed to change. Hence, one finds 
that the pressure is a contribution of two polytropes: a n = 3 polytrope for the mixture of 
gas and radiation and a n = 3/2 polytrope for the NC contribution, which is considered as a 
perturbation. The question of stability can be understood as follows. Defining V = 1 + -, a 
n = 3 (r = 4/3) polytrope represents a marginally stable configuration [33]. If one introduces 
corrections due to general relativity, one shows that, in fact, this polytrope is unstable. One 
follows here the purely Newtonian analysis [33]. The total energy is written as the sum of the 
internal energy (~ PV) and the gravitational potential energy (~ ^ ). The pressure is given 
by P = K\p Tl + XK2p^ 2 and to ensure that the second term is a perturbation A-^p r2 ~ ri <C 1. 
The energy can be written as 



E = k PV - h 



GM 2 
R 



CxMpl 1 - 1 + C 2 Mp 



r 2 -i 



hM^p\l\ 



(31) 



where ko, kx , k$ , C% , C 2 are constants and p c is the central density. The value of M at the critical 
point (§jf) = is computed and one performs 



dlnM 



2 V 3 



A 



3C 2 (r 2 - 1) 



(r 2 -ri) P ; 



r 2 -ri 



dlnp c 2V 1 3 J ' ' 2Ci(Ti - 1) 
A necessary, but not sufficient, condition for stability of the configuration is that 



(32) 



>0, 



for M given by If = [35]. Here, Y x = 4/3 and T 2 = 5/3 > Ti, so > °- Thus > the effect 



of noncommutativity is to turn this configuration into a stable one. 



To quantify the relevance of noncommutativity for these type of stars one chooses the Sun, 
as it is a very typical main-sequence star. In the model investigated here, at the center, 
1 - p s « 10" 3 , p c = 1.53 x 10 5 kg/m 3 and p N = 0.829 [Sj- The ratio between the 
noncommutative and the commutative contributions to the pressure is 

PNC \K 2 p 5/3 QfiRvln -16 x { oo\ 

~Pc~ = ~K^yW = 9M X 10 A ' (33) 

where Pnc = ^K 2 p 5 ^ 3 is the noncommutative correction to the pressure, Pc = K\p^^ is the 
usual (commutative) term and A must be expressed in S.I. units. One computed this value at 
the center of the star using Eqs. (|29|) and ([30]) . The value of A will be discussed later (Section 
[5|), but one sees that this ratio is fairly small. 

4-2. White Dwarfs 

As one has seen, for main-sequence stars the noncommutative correction is rather small and, even 
though it plays a role in the stability of the star, one should seek more energetic configurations. 
White dwarfs are stars that have masses comparable to th Sun and planet sizes [39|. At 
their center, the range of possible mass densities is 10 8 kg/m 3 < p c < 10 12 kg/m 3 . At these 
densities the matter is completely ionized and the electrons are free. The treatment of matter 
must necessarily involve quantum mechanics. The simplest model of white dwarfs consists 
in considering an ideal electronic gas maintained in equilibrium by gravity. Although the 
central temperature is T ~ 10 K, it is much smaller than the temperature associated with 
the Fermi energy (T ~ 10 10 K). This justifies describing white dwarfs, at first approximation, 
as a degenerate electron gas. 

Defining the electronic molecular weight as fx e = p/n e m p , where n e is the electronic particle 
number density and m p is the proton mass. One can rewrite Eq. (|25p as 

, = lt ^±^ = (*) V3 ±- (-£-] = !0-= (JL) V3 , (34) 
m e c \m p J m e c\p e J \p e J 

where m e is the mass of the electron and this result is presented in S.I. units. One considers 
p e = 2, as this corresponds to a star composed mainly by Helium. For the possible values of 
the central density, one chooses a low density value p = 10 8 kg/m 3 , an intermediate density 
(pf ~ w e c) value p = 2 x 10 9 kg/m 3 and a high one p = 10 12 kg/m 3 . One finds xg = 0.37, 
xg = 1.01, x±2 = 7.99, respectively. One now quantifies the relevance of noncommutativity, 
using Eq. (f2"7|) : 

Pnc _ 48x 5 2 

Pc "l5/(x) AmeC - (35) 

This function is plotted in Fig. [2] for the usual values of x found in white dwarfs. One 
sees that the relevance of noncommutativity is just a few Am e c 2 . Since m e c 2 ~ 0.5 MeV and 
A ~ {Eqc)~ 1 (see Section [5]), this correction is much smaller than other improvements to the 
simplest white dwarf model, such as Coulomb interaction for low densities, and general relativity 
for high densities. Also neglected are the corrections due to rotation and magnetic fields, which 
are likely to be greater than the noncommutative correction. 

Let us now analyse the problem of stability (for details see e.g. Ref. [34J ) : the 
noncommutative correction is proportional to x 5 tx /5 5 / 3 (Eq. (|34p ) which represents a V = 5/3 
polytrope. This lies in the region of stability and hence noncommutativity does not bring any 
instability for white dwarfs. 




Figure 2. Ratio between the noncommutative and the commutative contributions to the 
pressure in units of Am e c 2 for white dwarfs and neutron stars. 



4-3. Neutron Stars 

Consider now that p 10 12 kg/m 3 , a density for which the protons and electrons combine 
and through the weak interaction form neutrons. These objects must be described by general 
relativity since ~ 0.1. Their masses are about the solar mass (M ~ 1-M©) and their radius 
about R ~ 10 km. 

These stars are mainly composed by neutrons, hence the simplest model for a neutron star, 
the Oppenheimer-Volkoff (OV) model [ID] , considers an ideal neutron gas counterbalanced by 
gravity. This model is rather unrealistic as one must introduce nuclear interaction effects. 
However, this simpler formulation has some advantages, the first being that one can compute 
the thermodynamical quantities in a closed and analytical form. Furthermore, the obtained 
results do not differ significantly from a more realistic EoS [38J; finally the OV EoS is stiffer 
than other EoS that include nuclear interaction [51]. This neutron gas is degenerate since the 
Fermi energy of neutrons (~ 1 GeV) is much greater than the thermal energy associated to the 
gas (~ 10" 1 MeV). 

Using Eqs. ([25]), ([26]) and ([27]), one obtains 
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- = ( - ^ + (Am - c2) yJ ' (37) 

P = ^7^(^ + (A-nC 2 )^, (38) 

where m n is the neutron mass. One cannot obtain the EoS P = P(u) in an analytical form. 
The relevance of noncommutativity is depicted in Fig. [2] with the modification m e — > m n . For 
the usual baryonic matter, n = 0.15 x 10 45 m~ 3 [38J, one obtains x = 0.35 and inserting this 
value into Eq. (|38p . one obtains that -^pg 2 = 2.1(Am n c 2 ). Notice that as m n ~ 10 3 m e , the effect 
of noncommutativity is much more important for this class of stars than for white dwarfs. 

For neutron stars the question of stability is more complex since general relativity is required. 
It involves solving numerically the Oppenheimer-Volkoff equation with the OV EoS. 



5. Discussions and conclusions 

In this contribution, the most relevant effects of noncommutativity for astrophysical objects 
were estimated. Noncommutativity is introduced through a deformed dispersion relation and 



the relevant thermodynamical quantities were calculated using the grand canonical ensemble. 
These results are applied to physical models describing main-sequence stars, white-dwarfs and 
neutron stars. 

In what concerns the noncommutative parameter, A, (cf. Eqs. ([1]) and Q), an upper bound 
for its value can be inferred from the breaking of Lorentz invariance applied to the problem of 
ultra-high energetic cosmic rays. The most stringent limit is A < 2.5 x 10 -19 GeV -1 = 1.6 x 1CP 9 
J -1 . For the stars, one is certainly in the low-energy limit, and hence this justifies the expansion 
of the quantities up to first order in A. 

Using A ~ 10 -19 GeV -1 , one can estimate the ratio between Pnc an d Pc f° r the main- 
sequence stars (as an example one chooses the Sun), white-dwarfs and neutron stars. One 
obtains ^ ~ 10" 25 , ~ 1(T 22 and tyg- ~ 1(T 19 (cf. Eqs. ([Ml), QSE) and pgjlh respectively. 
These results indicate that noncommutative correction is fairly small, actually smaller than 
other usually neglected standard effects in stellar physics. As one expects the relevance of 
noncommutativity to grow as one considers denser configurations, one is lead to conclude that 
noncommutative might be relevant for black holes. In fact, full phase-space noncommutativity 
is shown to be quite important for Schwarzschild black holes [42]. 

For main sequence stars and white-dwarfs, the question of stability has been analysed and 
shown that noncommutativity does not introduce any instability. Actually, for main sequence 
stars, noncommutativity is beneficial for the star stability. 
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